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1. Introduction 

Kinetic theory is an important part of general statistical physics related 
to phenomena which cannot be understood on the thermodynamic or fluid 
models level. First of all wc mean (local) fluctuations from equilibrium state 
and a lot of complex phenomena ^. Also, it is well-known in plasma physics, 
e.g., that only a kinetic approach can describe Landau damping, intra-beam 
scattering, while Schottky noise and associated cooling techniques require 
understanding of the spectrum of local fluctuations of the beam charge 
density ^. In this paper we consider the applications of a new numeri- 
cal/analytical technique based on wavelet analysis approach for calculations 
related to the description of complex behaviour in the framework of the 
general BBGKY-hierarchy ^. We restrict ourselves to rational/polynomial 
type of nonlinearities (with respect to the set of all dynamical variables 
(n-particle distribution functions in the case under consideration)) which 
allows us to use our results from Refs. which are based on the appli- 
cation of wavelet analysis technique and variational formulation of initial 
nonlinear (pseudodifferential) problems. 

Wavelet analysis is a set of mathematical methods which give us a pos- 
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sibility to work with well-localized bases in 

functional spaces and provide maximum sparse forms for the general 
type of operators (differential, integral, pseudodifferential) in such bases. 
It provides the best possible rates of convergence and minimal complexity 
of algorithms inside and, as a result, saves CPU time and HDD space 
Our main goals are an attempt of classification and construction of 
possible nontrivial states in the system imder consideration. First of all 
we are interested in the following states: localized, chaotic-like patterns, 
localized (stable) patterns. We start from the corresponding definitions (at 
this stage these definitions have only qualitative character). 

1. By localized state (localized mode) we mean the corresponding (par- 
ticular) solution of the system under consideration which is localized in 
maximally small region of the phase space. 

2. By chaotic pattern we mean some solution (or asymptotics of solu- 
tion) of the system under consideration which has equidistribution of energy 
spectrum in a full domain of definition. 

3. By localized pattern (waveleton) we mean (asymptotically) stable 
solution localized in relatively small region of the whole phase space (or a 
domain of definition). In this case all energy is distributed during some 
time (sufficiently large) between few localized modes (from point 1) only. 
We believe, it is a good model for a plasma fusion (energy confinement) in 
the nearest future. It is also obvious that such states are very important in 
many areas of statistical physics. 

In all cases above, by the system imder consideration we mean the full 
BBGKY hierarchy or some cut-off of it. Our construction of cut-off of the 
infinite system of equations is based on some criterion of convergence of 
the full solution. This criterion is based on a natural norm in a proper 
functional space, which takes into account (non-perturbatively) the under- 
lying multiscale structure of complex statistical dynamics. According to 
our approach the choice of the underlying functional space is important 
for understanding the corresponding complex dynamics. It is obvious that 
we need to fix accurately the space in which we construct the solutions, 
evaluate convergence etc. and in which the very complicated infinite set 
of operators acts which appears in the BBGKY formulation. We underline 
that many concrete features of the complex dynamics are related not only 
to the concrete form/class of operators/equations but depend also on the 
proper choice of function spaces, where operators act ^. It should be noted 
that the class of smoothness of the functions under consideration plays a 
key role in the following (more details will be considered elsewhere). 

In Sec. 2 the kinetic BBGKY-hierarchy is formulated and an important 
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particular case is described. In Sec. 3 we present the explicit analytical 
construction of solutions of the hierarchy, which is based on tensor alge- 
bra extensions of bases generated by the hidden multiresolution structure 
and proper variational formulation leading to an algebraic parametrization 
of the solutions. We give the explicit representation of the hierarchy of 
N-particle reduced distribution functions in the basis of highly-localized 
generalized coherent (regarding the underlying afhnc group) states given 
by polynomial tensor algebra of wavelets, which takes into account contri- 
butions from all underlying hidden multiscales, from the coarsest scale of 
resolution to the finest one, to provide full information about stochastic 
dynamical process. So, our approach resembles Bogolyubov's and related 
approaches but we don't use any perturbation technique (like virial expan- 
sion) or linearization procedures. Numerical modeling as in general case 
as in particular cases of the Vlasov-likc equations shows the creation of 
various internal (coherent) structures from localized modes, which are re- 
lated to stable (equilibrium) or unstable/chaotic type of behaviour and the 
corresponding pattern (waveletons) formation. 



2. BBGKY-Hierarchy 

Let M be the phase space of an ensemble of N particles (dimM = 67V) 
with coordinates Xi = {qt,Pi), i = 1, N, qi = {ql,q^, ) e i?^, Pi = 
{pjjPijPi) ^ 9 = ■ • • € R^^ ■ Individual and collective mea- 

sures are: /Zj = dxi — dqidpi, fi = Ili^iMi- 'I'^e distribution fimction 
Dn{xi, . . . ,XN;t) satisfies Liouville's equation of motion and the normal- 
ization constraint for an ensemble with the Hamiltonian ■ 

^^^ = {Hn,Dn} J Dn{xi, . . . ,XN;t)dfi = I. (1) 

Our constructions can be applied to the following general Hamiltonians: 

N 2 

^^ = E(^ + f^'(9))+ E Uij{qi,qj), (2) 

i=l l<i<j<N 

where the potentials Ui{q) = Ui{qi, . . . , qN) and Uij{qi, qj) are restricted to 
rational functions of the coordinates. 

Let Ls and Lij be the Liouvillean operators (vector fields) 

= E (1^ - - E ^^^^ (3) 

''^ dqi dpi dqj dpj' 
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Let F/v(a;i, . . . ,XN;t) = J2sn Dn{xi, . . . ,XN;t) he the A'^-particle distribu- 
tion function {Sn is permutation group of N elements). Then we have the 
hierarchy of reduced distribution functions {V is the volume) 



Fs{xi, ...,Xs;t) = V' 



j DN{xi,...,XN]t] l^i (5) 

^11 ^ \r 



s-V\<i<N 

After standard manipulations we arrive at the BBGKY-hierarchy ^: 



dt 



+ LsFs = ^ j d/Us+i (6) 



It should be noted that we may apply our approach even to more general 
formulation (nonlinear) than (6). For s=l,2 we have, from the general 
BBGKY-hierarchy: 



dFi{xi;t) _^ P]__d_ 

dt m dqi ' V 



Fi{xi;t) = ^ j Ax2Li2F2{xi,X2;t), 
dF2{x\,X2\t) ^ (Pi d P2 d 



dt 



{'-^ + '-^-Li.)F2{xi,X2;t) (7) 
V m oqi m oq2 ) 



= yj dx3{Li3 + L23)F3{xi,X2;x3;t). 



Ill most cases, one is interested in a representation of the form 
Fk{xi,...,Xk;t) = Yii^i^iixi't) + Gk{xi,...,Xk;t), where Gk are cor- 
relators. Additional reductions often lead to simplifications, the simplest 
one, Gk = 0, corresponding to the Vlasov approximation. Such physically 
motivated ansatzes for Fk formally replace the linear (in Fk) and pseudod- 
ifferential (in general case) infinite system (6), (7) by a finite-dimensional 
but nonlinear system with polynomial nonlincaritics (more exactly, miilti- 
linearities) Our key point in the following consideration is the proper 
generalization of the perturbative multiscale approach of Bogolyubov. 



3. Multiscale Analysis 

The infinite hierarchy of distribution functions satisfying system (6) in the 
thermodynamical limit is: 

F = {Fq, Fi{xi;t), F2{xi,X2; t),..., Fn{xi, . . . ,XN;t), . . .}, (8) 

where Fp{xi, . . . ,Xp;t) e RP, H° = R, RP = L^{R^p) (or any different 
proper functional space), F e H°° = H° ® H'^ ® . . . ® HP ® . . . with 
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the natural Fock space like norm (guaranteeing the positivity of the full 
measure): 

{F, F) = F^ + Y,f F^ixi, •.•,Xi;t)f[ M^ (9) 

First of all we consider F = F{t) as a function of time only, F G L'^{R), via 
multiresolution decomposition which naturally and efficiently introduces 
the infinite sequence of the underlying hidden scales Because the affine 
group of translations and dilations generates multiresolution approach, this 
method resembles the action of a microscope. We have the contribution 
to the final result from each scale of resolution from the whole infinite 
scale of spaces. Wc consider a multiresolution decomposition of i^(-R) ^'^ 
(of course, we may consider any different and proper for some particular 
case functional space) which is a sequence of increasing closed subspaces 
Vj e L'^{R) (subspaces for modes with fixed dilation value): 

...F_2 C V-i C Vb C C F2 C ... (10) 

The closed subspace Vj{j e Z) corresponds to the level j of resolution, or to 

the scale j and satisfies the following properties: let Wj be the orthonormal 
complement of Vj with respect to Vj+\ : Vj+\ = Wj . Then we have 

the following decomposition: 

oo 

{F{t)}= Wj or {F{t)} = Vo^Wj, (11) 

— oo<j<oo i— 

in case when Vq is the coarsest scale of resolution. The subgroup of trans- 
lations generates a basis for the fixed scale number: span^^^l^^'^^^l^''^ ~ 
fc)} = Wj. The whole basis is generated by action of the full affine group: 

spanfce^,,.g^{2^/2^'(2^t - k)} = span^^^.^^f^^.a = {F{t)}. (12) 

Let the sequence {Vj},Vj C L^{R) correspond to multiresolution analysis 
on the time axis, {V^^'} correspond to multiresolution analysis for coordi- 
nate Xi, then V"^^ = V^^ ® . . . ® V^^" ® Vj corresponds to the multires- 
olution analysis for the n-particle distribution function Fn{xi, . . . ,Xn',t). 
E.g., for n = 2: = {/ : fixi,X2) = Y^kuk2°-kiM4>^{xx - ki,X2 - 

ki), ak,M £ ^'^(Z^)}, where (p{xi,x2) = (/)^(xi)</>2(a;2) = ^ (j)'^{xi,X2), 
and 4>^(xi) = (t>{xi) form a multiresolution basis corresponding to {Vj'^}. 
If {(^^(xi — £)}, £ & Z form an orthonormal set, then (p'^{xi — ki,X2 ~ fc2) 
form an orthonormal basis for Vq. So, the action of the affine group gener- 
ates multiresolution representation of LF'{E?'). After introducing the detail 
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spaces Wj, we have, e.g. = Vq © Wq. Then the 3-coinponent basis for 
Wq is generated by the translations of three functions 

Also, we may use the rectangle lattice of scales and one-dimensional wavelet 
decomposition: 

f{xi,X2) = J2 *i,fe)*i,^® *i,fe(a;i,a;2), 

where the basis functions ^i^e (8) ^j^k depend on two scales 2~* and . 
After constructing the multidimensional basis wc may apply one of the 
variational procedures from Refs. We obtain our multiscale/multire- 

solution representations (formulae (17) below) via the variational wavelet 
approach for the following formal representation of the BBGKY system 
(6) (or its finite-dimensional nonlinear approximation for the n-particlc 
distribution functions) with the corresponding obvious constraints on the 
distribution functions. 

Let L be an arbitrary (non)linear differential/integral operator with 
matrix dimension d (finite or infinite), which acts on some set of func- 
tions from L'^{Q'^"): * = "^{t, xi, X2, ■ ■ ■) = (^*^(t, Xi, X2, ...),.. ., 

"^'^{t, Xi,X2, ■ ■ ■)), Xi G ^ C R®, n is the number of particles: 



L^ = L{Q,t,Xi)^{t,Xi) = 0, (13) 

where 

Q = <3do,di,d2,...(*i xi,X2,..., d/dt, d/dxi, d/dx2, j Mfc) = 

''"'V^'"" / s(9y»( d Nn/ d f ,,,, 

_1^onn...it,XuX2,...)[g-^) (^j (^j ...yMfe.(14) 



ia,ii,i2.... = l 



Let us consider now the N mode approximation for the solution as the 
following ansatz: 



JV 



*^(i,a;i,a;2,...) = ^ aigi,i^„Ai^<S,Bi^iSiCi^...{t,xi,X2,...).{15) 



=1 



We shall determine the expansion coefHcients from the following conditions 
(different related variational approaches are considered in ^~^): 



lk,M,- = J {L^i>^)Ako{t)BkAxi)Ck,{x2)dtdx,dx2... = 0. 



(16) 
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Thus, we have exactly dN"' algebraical equations for diV" unknowns 
ciioM,---- This variational approach reduces the initial problem to the prob- 
lem of solution of functional equations at the first stage and some algebraical 
problems at the second. We consider the multircsolution expansion as the 
second main part of our construction. So, the solution is parametrized 
by the solutions of two sets of reduced algebraical problems, one is linear 
or nonlinear (depending on the structure of the operator L) and the rest 
are linear problems related to the computation of the coefficients of the 
algebraic equations (16). These coefficients can be found by some wavelet 
methods by using the compactly supported wavelet basis functions for the 
expansions (15). As a result the solution of the equations (6) has the follow- 
ing multiscale or multircsolution decomposition via nonlinear high-localized 
eigenmodes 

F{t,Xi,X2,...) = ^ ttijW lSiV^{t,Xi,X2,...), 

Vi{t)=V'/''''"{t) + Y,Vi\wit), coir. 2', (17) 

1>N 

U\X,) = U'^'"^{Xs) + ULik^n^s), k^m ^ 2", 

m>M 

which corresponds to the full multircsolution expansion in all underlying 

time/space scales. The formulae (17) give the expansion into a slow part 
*?^°M ^"^d f^st oscillating parts for arbitrary A'^, M. So, we may move from 
the coarse scales of resolution to the finest ones for obtaining more detailed 
information about the dynamical process. In this way one obtains contri- 
butions to the full solution from each scale of resolution or each time/space 
scale or from each nonlinear eigenmode. It should be noted that such repre- 
sentations give the best possible localization properties in the correspond- 
ing (phase)space/time coordinates. Formulae (17) do not use perturbation 
techniques or linearization procedures. Numerical calculations are based 
on compactly supported wavelets and related wavelet families and on 
evaluation of the accuracy on the level N of the corresponding cut-off of 
the full system (6) regarding norm (9): 

\\F^+^ -F^\\<s. (18) 

To summarize, the key points are: 

1 . The ansatz-oriented choice of the (multidimensional) basis related to 
some polynomial tensor algebra. 

2. The choice of a proper variational principle. A few projection (Galer- 
kin-like) principles for constructing (weak) solutions are considered in 
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The advantages of formulations related to biorthogonal (wavelet) decom- 
positfon should be noted. 

3. The choice of basis functions in the scale spaces Wj. They correspond 
to highly-localized (nonlinear) oscillations/excitations, nontrivial local (sta- 
ble) distributions/fluctuations, etc. Besides fast convergence properties we 
note the minimal complexity of all underlying calculations, especially by 
choosing wavelet pac;kc!ts which minimize Shann(m's entropy. 

4. Operator representations providing maximum sparse representations 
for arbitrary (pseudo) differential/ integral operators d//da;, d"//da;", 
jTix,y)f{y)dy),etc. 

5. (Multi)linearization. Besides the variation approach we can con- 
sider also a different method to deal with (polynomial) nonlinearities: para- 
products-like decompositions 

4. Example: Vlasov equation 

As a particular case we consider the Vlasov approximations F2 = FiFi, 
F3{xi,X2,x:i) = Pi{xi)F2{xj,Xk) in Eqs. (7), which are important in 
plasma physics. 

This is a particular case of the general form (13) (for simplicity we 
consider only one variable) 

Q{R, x)^{x) = P{R, x)^{x) or = L{R, x)^{x) = 0, (19) 

where R = R{x,d/dx,^) is not more than a rational (operator) function. 
We have the following representation for the N mode approximation for 
the solution of the Vlasov equation via expansion in some high-localized 
wavelet-like basis (other independent variables are considered analogously): 

N 

^''{x)=J2<M^)- (20) 

We shall determine the expansion coefficients from the following variational 
conditions: 

L^ = J{L^^)M^)dx = 0. (21) 

We have exactly dN algebraical equations for dN unknowns a^. So, 
the variational approach reduces the initial problem (7) to the problem of 
the solution of functional equations at the first stage and some algebraical 
problems at the second stage. As a result we have the following reduced 
algebraical system of equations (RSAE) on the set of unknown coefficients 
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ooortinate (q) " « momentum (p) 



Figure 1. Localized mode contribution to distribution function. 




coortinate(q) " " momentum (p) 



Figure 2. Chaotic-like pattern. 



af^ of the expansion (20): 

H{Qij,af,ai) = M{Pij,a^ , pj), (22) 

where the operators H and M are algebraizations of the RHS and LHS 

of the initial problem (19). Qij (Pij) are the coefficients of LHS (RHS) 
of the initial system of differential equations (19) and, as consequence, 
are the coefHcients of the RSAE. I = (ii, ig_|_2), J = {ji, ■■■,jp+i) are 
multiindexes labeling aj and /?/, the other coefficients of (22) are 

Pj = {Pn-j.^A^ [ n (23) 
i<jfc<P+i 

where p is the degree of the polynomial operator P, 

il,---,iq+2 
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and q is the degree of the polynomial operator Q, = (1, q + 2), t^j^ = 

We may extend our approach to the case when there are additional con- 
straints on the set of dynamical variables and additional averaged terms 
also. In these cases, by using the method of Lagrangian multipliers one 
again may apply the same approach, but for the extended set of variables. 
As a result, one obtains the expanded system of algebraical equations, anal- 
ogous to the system (22). Then, one again can extract from its solution the 
coeflacients of expansion (20). It should be noted that if one considers only 
the truncated expansion (22) with N terms one has a system of N xd alge- 
braical equations of degree i = max{p, q} which coincides with the degree 
of the initial system. 




coordinate (q) ° » momentum (p) 

Figure 3. Localized waveleton pattern. 



5. Demonstration of typical patterns 

Formulae (15), (17), (20) provide, in principle, a fast convergent decomposi- 
tion for the general solutions of the systems (6), (7) in terms of contributions 
from all underlying hidden internal scales. Of course, we cannot guarantee 
that each concrete system (6), (7) with fixed coefficients will have a priori 
a specific type of behaviour, either localized or chaotic. Instead, we can 
analyze if typical structures described by definitions 1-3 are present. To 
classiiy the qualitative behaviour we apply standard methods from general 
control theory or really use the control ^. We will start from a priori un- 
known coefficients, the exact values of which will subsequently be recovered. 
Roughly speaking, we will fix only class of nonlinearity (polynomial in our 
case) which covers a broad variety of examples of systems (6), (7). 



12 



As a simple model we choose band-triangular non-sparse matrices (ay) 
from (17) or (20) in particular cases. These matrices provide tensor struc- 
ture of bases in (extended) phase space and are generated by the roots of 
the reduced variational (Galerkin-like) systems (16) or (21), (22). As a sec- 
ond step we need to restore the coefficients of (6), (7) from these matrices 
by which we may classify the types of behaviour. We start with the local- 
ized mode, which is an "elementary" eigenfunction, Fig. 1, corresponding 
to def. 1, which was constructed as a tensor product of the two Daubechies 
functions. Fig. 2, corresponding to def. 2, presents the result of summation 
of series (17) up to value of the dilation/scale parameter equal to six on the 
bases of symmlets with the corresponding matrix elements equal to one. 
The size of matrix is 512x512 and as a result we provide modeling for one- 
particle distribution function corresponding to standard Vlasov like-system 
(7) with F2 = Ff. So, different possible distributions of the root values of 
the generical algebraical system (22) provide qualitatively different types 
of behaviour. The above choice provides us by chaotic-like equidistribution 
distribution. But, if we consider a band-like structure of matrix (a^) with 
the band along the main diagonal with finite size (<C 512) and values, e.g. 
five, while the other values are equal to one, we obtain localization in a 
fixed finite area of the full phase space, i.e. almost all energy of the system 
is concentrated in this small volume. This corresponds to definition 3 and is 
shown in Fig. 3, constructed by means of DaTibechic^s-based wavelet pack- 
ets. Depending on the type of solution, such localization may be present 
during the whole time evolution (asymptotically-stable) or up to the needed 
value from time scale (e.g. enough for plasma fusion/confinement). 

Now we discuss how to solve the inverse/synthesis problem or how to 
restore the coefficients of the initial systems (6), (7). Let 



be the system (13) with the fixed coefficients The corresponding 

solution is represented by formulae (15), (17) or (20), which are 
parametrized by roots of reduced algebraic system (22) and constructed 
by some choice of the tensor product bases from part 1. The proper coun- 
terpart of the system (25) with prescribed behaviour corresponding 
to a given choice of both tensor product structure and coefficients {fly } 
described above, corresponds to the class of systems like (13) but with 
undetermined coefficients and has the same form 



Our goal is to restore coefficients from (25), (26) and explicit represen- 
tations for solutions and This is a standard problem in the adaptive 



L0(Q°)*° = 



(25) 



L"(Q")*" = 0. 



(26) 
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control theory one adds a controlling signal u{x,t) which deforms the 
controlled signal ^{x, t) from the fixed state \E'°(x, t) to the prescribed one 
^'"(a;, t). At the same time one can determine the parameters Q" ^ ■ Finally, 
we apply two variational constructions. The first one gives the systems of 
algebraic equations for unknown coefficients, generated by the following set 
of functionals 

j ({1° - L^)^l, K)dl^N, (27) 

where A'' means the A-order approximation according to formulae (15). 

The unknown parameters Q* are given by Q* = limAr^oo Q]v- '^^^ second 
is an important additional constraint on the region /xq in the phase space 
where we are interested in localization of almost all energy E = J (\I/")d/Lt, 
where E is the proper energy functional ^ (Marsden-like). 

We believe that the appearance of nontrivial localized patterns observed 
by these simple methods is a general effect which is also present in the full 
BBGKY-hierarchy, due to its complicated intrinsic multiscale dynamics 
and it depends on neither the cut-off level nor the phenomenological-like 
hypothesis on correlators. So, representations like (17) and the prediction of 
the existence of the (asymptotically) stable localized patterns/states (en- 
ergy confinement states) in BBGKY-like systems are the main results of 
this paper. 
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